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Abstract: We use recently proposed method of ratios to assess the quality of geometrical 
scaling in deep inelastic scattering for different forms of the saturation scale. We consider 
original form of geometrical scaling (motivated by the Balitski-Kovchegov (BK) equation 
with fixed coupling) studied in more detail in our previous paper, and four new hypotheses: 
phenomenologically motivated case with dependent exponent A that governs small x 
dependence of the saturation scale, two versions of scaling {running coupling 1 and 2) that 
follow from the BK equation with running coupling, and diffusive scaling suggested by 
the QCD evolution equation beyond mean field approximation. It turns out that more 
sophisticated scenarios: running coupling scaling and diffusive scaling are disfavored by 
the combined HERA data on e'^p deep inelastic structure function F2. 
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1 Introduction 

Geometrical scaling (GS) has been introduced in Ref. [1] in the context of low x Deep 
Inelastic Scattering (DIS). It has been conjectured that 7*p cross-section cr^*p(a;, Q^) = 
47r^ 06111-^2(3; 5 Q'^)/Q'^ which in principle depends on two independent kinematical variables 
and W (i.e. 7*p scattering energy), depends only on a specific combination of them, 
namely upon 

called scaling variable. Bjorken x variable is defined as 



(1.2) 



and Mp denotes the proton mass. In Ref. [1], following Golec-Biernat-Wiisthoff (GBW) 
model [2], function Qs{x) - called saturation scale - was taken in the following form 

Qli^) = Ql{^^ ' ■ (1-3) 

Here Qq and xq are free parameters which can be extracted from the data within some 
specific model of DIS, and exponent A is a dynamical quantity of the order of A ~ 0.3. In 
the GBW model Qq = I GeV/c and = 3 x IQ-l 

In our previous paper [3] (see also [4]) we have proposed a simple method of ratios 
to assess in the model independent way the quality and the range of applicability of GS 
for the saturation scale defined in Eq. (1.3). Here we follow the same steps to test four 
different forms of the saturation scale that have been proposed in the literature. 

Geometrical scaling is theoretically motivated by the gluon saturation phenomenon 
(for review see Refs. [5, 6]) in which low x gluons of given transverse size ~ l/Q^ start 
to overlap and their number is no longer growing once is decreased. This phenomenon 
- called gluon saturation - appears formally due to the nonlinearities of parton evolution 
at small x given by so called JIMWLK hierarchy equations [7] which in the large A^c 
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limit reduce to the Balitsky-Kovchegov equation [8] . These equations admit travehng wave 
solutions which explicitly exhibit GS [9]. An effective theory describing small x regime is 
Color Glass Condensate [10, 11]. 

Gluon saturation takes place for Bjorken x much smaller than 1. Yet in Ref. [3] we have 
shown that GS with saturation scale defined by Eq. (1.3) works very well up to much higher 
values of x, namely up to a; ~ 0.1. In this region GS cannot be attributed to the saturation 
physics alone. Indeed, it is known that GS scaling extends well above the saturation scale 
both in the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi [12] (DGLAP) [13] and Balitsky- 
Lipatov-Fadin-Kuraev [14] (BFKL) [15] evolution schemes once the boundary conditions 
satisfy GS to start with. It has been also shown that in DGLAP scheme GS builds up 
during evolution for generic boundary conditions [16] . Therefore in the kinematical region 
far from the saturation regime where, however, no other scales exist {e.g. for nearly massless 
particles) it is still the saturation scale which governs the behavior of the j*p cross-section. 

The form of saturation scale given by Eq.(1.3) is dictated by the asymptotic behavior 
[9] of the Balitsky-Kovchegov (BK) equation [8], which is essentially the BFKL equation 
[14] supplied with a nonlinear damping term. It has been first used in the papers by K. 
Golec-Biernat and M. Wiisthoff [2] where the saturation model of inclusive and diffr active 
DIS has been formulated and tested phenomenologically. 

Since the original discovery of GS in 2001 there have been many theoretical attempts to 
find a "better" scaling variable which is both theoretically justified and phenomenologically 
acceptable. An immediate generalization of the saturation model of Refs. [2] has been done 
in Ref. [17] where DGLAP [12] evolution in has been included. Although the exact 
formulation of DGLAP improved saturation model requires numerical solution of DGLAP 
equations, one can take this into account phenomenologically by allowing for an effective 

dependence of the exponent A = Aphn(Q^) which is indeed seen experimentally in the 
low x behavior of F2 structure function (see e.g. Refs. [17, 18] and Fig. 1). This piece 
of data can be relatively well described by the linear dependence of Aphn(Q^) on logQ^ 
leading to the scaling variable of the following form 

rphn = QV«+^'°^«'/«' (1.4) 

In another approach to DIS at low x one considers modifications of BK equation 
through an inclusion of the running coupling constant effects. Depending on the approxi- 
mations used two different forms of scaling variable have been discussed in the literature 
[9]: 

r,,i = g2g-^.Viog(iA) (1.5) 

and [19] 

rrc2 = gV/i°g(«VQ2) (1.6) 

where subscripts " rc" refer to " running coupling" . Note that from phenomenological point 
of view (1.6) is in fact a variation of (1.4) where a different form of dependence has 
been used. Finally, generalization of the BK equation beyond a mean-field approximation 
leads to so called diffusive scaling [20] characterized by yet another scaling variable: 
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These different forms of scaling variable (except (1.4)) have been tested in a series of 
papers [21-23] where the so called Quality Factor (QF) has been defined and used as a 
tool to assess the quality of geometrical scaling. In the following we shall use the method 
developed in Refs. [3, 4] to test hypothesis of GS in scaling variables (1.4)-(1.7) and to 
study the region of its applicability using combined analysis of e~^p HERA data [24] . We 
shall also compare our results with earlier findings of Refs. [21-23]. 

Our results can be summarized as follows: more sophisticated scenarios i.e. running 
coupling scaling and diffusive scaling are disfavored by the combined HERA data on e~^p 
deep inelastic structure function F2. In contrast, phenomenologically motivated case with 
dependent exponent A and the originally proposed form of the saturation scale [1] with 
fixed A exhibit high quality geometrical scaling over the large region of Bjorken x up to 
0.1. The fact that GS is valid up to much larger Bjorken x's than originally anticipated 
has been already used in an analysis of GS in the multiplicity px spectra in pp collisions 
[25]. 

In Sect. 2 we briefly recapitulate the method of ratios of Ref. [3] and define the criteria 
for GS to hold. In Sect. 3 we present results for 4 different scaling variables introduced 
in Eqs.(1.4)-(1.7). Finally in Sect. 4 we compare these results with our previous paper [3] 
and with the results of Refs. [21-23]. 



2 Method of ratios 



Throughout this paper we shall use model-independent method used in Refs. [3, 4] which 
was developed in Refs. [26] to test GS in multiplicity distributions at the LHC. Geometrical 
scaling hypothesis means that 

Trpixi,Q^) = ^F{T) (2.1) 

where for simplicity we define 

(T'r*p[Xi,Q ) = ^2 • (^-^J 

Function F in Eq. (2.1) is a universal dimensionless function of r. In view of Eq. (2.1) 
cross-sections o''y*p{xi,Q'^) for different Xj's, evaluated not in terms of but in terms of 
T, should fall on one universal curve. This means in turn that if we calculate ratio of 
cross-sections for different Bjorken Xj's, each expressed in terms of r, we should get unity 
independently of r. This allows to determine parameter governing x dependence of r by 
minimizing deviations of these ratios from unity. Generically we denote this parameter as 
a, although for each scaling variable (1.4) - (1.7) it has a different meaning: a = f3, /j,, i' 
and K for Q^-dependent, running coupling (1 and 2) and diffusive scaling hypotheses, 
respectively. 

Following [3, 4] we apply here the following procedure. First we choose some Xref and 
consider all Bjorken Xj's smaller than Xy-^f that have at least two overlapping points in 
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(or more precisely in scaling variable r). Next we form the ratios 



aYp{xi,T{xi,Ql;a)) 
cr^*p(xi-ef5 T{xj-(,f, Q'j, j.g£; 0()) 



i?...w(a; r.) = , 02 (2-3) 



with 

Tfc = r(xi, Q^; a) = r(xrcf, Q^ref! «)• (2-4) 

By tuning a one can make Rxi.x^ci{(^',Tk) = 1 ± (5 for all r^. with accuracy of 5 for which 
following Ref. [3] we take 3%. 

For a ^ points of the same but different x's correspond generally to different 
r's. Therefore one has to interpolate the reference cross-section cr-y*p(xref, T(xref, Q^', a)) to 
Qk ref such that r(xref, Q^. j-gf, «) = Tfc as indicated in Eq. (2.4). This procedure is described 
in detail in Refs. [3, 4]. 

In order to find optimal value of parameter a that minimizes deviations of ratios (2.3) 
from unity we form the chi-square measure 

v2 (a) = - y{Rx.,..Ja;rk)-lf . . 

where the sum over k extends over all points of given Xi that have overlap with Xref and 
]^xi,x^gf is a number of such points. 

Finally, the errors entering formula (2.5) are calculated using 

ARx„x,Ja;Tkf = (2.6) 

'Aa^*p{Xi,T{Xi,Ql))\^ /Aa^*p(Xref,r(Xref,Q|ref))^ 



( ( n'^\\ I +( '~T T TW^^^SV ^ \ Rxi,x,^f{a;Tk)'^ + 6'^ 

where Aa^*p{T{x, Q^)) are experimental errors (or interpolated experimental errors) of 7*p 
cross-sections (2.2). For more detailed discussion of errors see Ref. [3]. 

In this way, for each pair of available Bjorken variables {xi, Xref), we compute the best 
value of parameter a, denoted in the following by a subscript min:^ amin(3^i) a^ref) and the 
corresponding x^- For GS to hold we should find a region in (xj, Xref) half-plane (note that 
by construction Xi < Xref) where amin(a^i) a^ref) is a constant independent of Xj and Xref, 
and the corresponding Xx^.x^cf small. 

We shall also look for possible violations of GS in a more quantitative way. In order 
to eliminate the dependence of aniin(a;, Xref) on the value of x, we introduce averages over 
X (denoted in the following by (...)) minimizing the following chi-square function: 

which gives the best value of (a) denoted as (aminC^^ref))- The sum in (2.7) extends over all 
x's such that amin(x,Xref) cxists and Nx^^i is the number of terms in (2.7). 



^Because it minimizes 
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Figure 1. Effective exponent Aphn from F2 at low x (3.2) from HERA and the linear fit of Eq. (3.3). 
Data points as in Ref. [17], see also [18]. 

Since GS is expected to work for small x's, the "average" value of scaling parameter 
(«mm(2^ref)) Supplies an information, up to what value of Xref GS is still working. For small 
Xref we expect (amin(2;ref)) to be constant, whereas for larger values we expect to see some 
dependence of (amin(iCref)) on Xref • A word of warning is here in order. Even if (amin(3;ref)) 
is a constant we have to look at the corresponding value of x^' too large obviously 
indicates violation of GS. 

To quantify further the hypothesis of geometrical scaling we form yet another chi- 
square function 



which we minimize to obtain ((amin(a^cut)))- 

Equation (2.8) allows us to see how well one can fit (amin(a;ref)) with a constant a up 
to Xref = a^cut- Were there any strong violations of GS above some xq, one should see a rise 
of ((amin(a;cut))) oncc Xcut becomes larger than xq. 

3 Results 

Let us now come back to the discussion of different scaling variables defined in Eqs. (1.4) - 
(1.7). All of them depend on one variational parameter, which we constrain analyzing 
ratios (2.3) for combined HERA e+p DIS data [24]. 

In the case of Q^-dependent exponent Aphn (1-4), however, there are in fact two pa- 
rameters, one of them (Aq) being fixed using our previous analysis of Ref. [3] where we 
have shown that GS scaling works very well with constant A = Aq: 




(2.8) 



Ao = 0.329 ±0.002. 



(3.1) 
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On the other hand looking at low x behavior of the F2 structure function it has been 
shown that [17, 18]: 

i^2(x,g2)~x-^phn(Q^) (3.2) 

where Aphn(Q^) can be well parametrized as 

Aphn(Q^) = 0.329 + 0.1 log(QV90) (3.3) 

(for in (GeV/c)^) as depicted in Fig. 1. Taking therefore scaling variable in the form of 
(1.4) with Ao = 0.329 we test in fact consistency of the slopes /3 as extracted from Fig. 1 
and by the procedure described in Sect. 2. Note that this is therefore a kind of perturbative 
two parameter fit, and as such it has a different status than the remaining Ansatze for the 
scaling variable (1.5) - (1.7). Similar remarks apply to the running coupling rc2 case (1.6), 
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Figure 3. Three dimensional plots of minimal values of functions (2.5) for different scaling 
variables: a) logarithmic scaling variable Tphn (1-4), running coupling b) rcl (1.5), c) rc2 (1.6) and 
d) diffusive scahng (1.7). 

where the scale of the logarithm has been fixed at 0.04 (following e.g. Ref. [22]). Then 
for all points > and t^c2 decreases with rising v. 

Let us first examine 3 dimensional plots of a^iriix, x^ci) (note again that a = {3, fi, u or 
AC, depending on the scaling variable). For GS to hold there should be a visible plateau of 
Omin over some relatively large part of (x, x^ef) space (recall that by construction x < Xref)- 
Looking at Fig. 2 one has to remember that the values of auim{x , x^ei) are subject to 
fluctuations that will be "averaged over" when we discuss more "integrated" quantities 
(ttmin) and ((amin))- Note that statistical errors of amm{x,Xrcf) which are quite large for 
small X are not displayed in Fig. 2. One can can conclude from Fig. 2 that for all 4 cases 
(1.4) - (1.7) there is rather strong dependence of amin{x, x^ef) for large values of x and x^ef- 
In the case of Q^-dependent scaling variable (1.4) (Fig. 2. a) and for the running coupling 
case (1.5,1.6) (Figs. 2.b, c), the values of parameters /3, fj, and v rise steeply for large x's. 
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Figure 4. Effective exponents (3.4) as functions of x for fixed — 10 GeV^/c^ (left) and as 
functions of for fixed x ~ 0.0001 (right). 



whereas for diffusive scaling parameter k is falling down rapidly. More closer look reveals 
that for running coupling rcl case (Fig. 2.b) there is in fact no distinct plateau, one can also 
see a systematic rise of /Xmin in a region of very small x's. Similarly for the diffusive scaling 
(Fig. 2.d) we see rather systematic growth of n^in for small x's with possible plateau in 
a small corner of very low x's. At first glance no plateau is neither seen for /3min(a^, a^ref) 
(Fig. 2. a). However ~ as will be shown in the following - because of considerable statistical 
uncertainties within the scale used in Fig. 2. a, very good description of GS with constant 
f3 is still possible. 

It is interesting to look at 3 dimensional plots of the corresponding values (2.5) 
shown in Fig. 3. Recal that for GS to hold one should observe small values of x^(Q^min) in 
the same region where amin is constant. This happens for Tphn (Fig. 3. a) where oscillates 
around 1 not exceeding 2 even for large values of x. Similarly Trci (Fig. 3.b) stays smaller 
than 2 up to X ~ 10~^ where jumps above 2. In this region, however, parameter /i is 
steadily decreasing with x. In contrast, in the case of Trc2 (Fig. 3.c) and Tds (Fig. 3.d) 
X^ have pronounced fluctuations and a plateau (if at all) is visible only below x ~ 10~^. 
However, in this region parameter v (corresponding to Fig. 3.c) rises with x, whereas k 
(corresponding to Fig. 3.d) exhibits rather strong fluctuations. 

Due to different functional dependence of the saturation scales entering Eqs. (1.4) - 
(1.7) variations of parameters f3, fj,, v and k differently influence pertinent scaling variable 
r. Therefore - before we turn to average quantities (. . .) and ((. . .)) displayed in Fig. 5 - 
let us deflne effective exponents Agg: 

Aeff(x,Q') = log(^)/log(x) (3.4) 

which depend on fitting parameters /3, fi, v and k. In Fig. 4 we plot these effective powers 
as functions of x and for the values of the parameters /3min, A*min, ^^min and fixed 
at the end of this Section. 

In order to find the scale relevant for a parameter entering definition of a given scaling 
variable r (1.4) - (1.7), for each scaling hypothesis separately we have varied this parameter 
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Figure 5. Averaged values (amin(a;rct)) (black squares) and ((amin(a;cut))) (red circles) for different 
scaling hypotheses: a) logarithmic effective exponent (1.4) with a — f3, running coupling scaling 
variables b) rcl (1.5) with a ~ fj, and d) rc2 (1.6) with a — v, and d) diffusive scaling (1.7) with 
a = K, respectively. 



around the best value by ±e and required that 

|Aeff(amin ± e;x,Q'^) - Acfr(amin; a::, Q^)| = 1 (3.5) 

for some typical values of x = 0.0001 and = 10 GeY^/c^. In this way in each case 
the value of e provides the reference scale for each variational parameter a = /3, /i, z/ or 
K. Therefore looking at Fig. 5 one should bear in mind that the span of the vertical axis 
corresponds to the variation of the effective exponent AAgff ~ ±1 around its best value. 

Looking at Figs. 5 we see immediately that the best scaling properties are exhibited 
by parameter (3 of Q^-dependent scaling variable Tphn (1.4). Parameter (3 is well described 
by a constant 

Po = ((/3min(0.08))) = 0.02 ± 0.001 (3.6) 

over 3 orders of magnitude in x. We have used the value of maximal Xcut = 0.08, since 
it was the value of Xcut for which Aq = 0.329 has been extracted in Ref. [3], although - 
as clearly seen from Fig. 5. a - GS in variable Tphn works well up to x ~ 0.2. There is an 
impressive agreement between both averages (/3min) and {{(3jam)), however the value (3.6) 
is five times smaller than expected from the fit to low x behavior of F2 structure function 
(3.3). 
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Figure 6. Averaged values (Aniin(a:rcf )) (black squares) and ((Aniin(a;cut))) (red circles) for scaling 
variable with constant exponent A (Eqs. (1.1) and (1-3)). Figure from Ref. [3]. 



For comparison in Fig. 6 we present the plot from Ref. [3] where (Amm) and ((Amin)) 
for scaling hypothesis with constant A (i.e. for /3 = 0) are shown. We see that the quality 
of a fit with a constant A is only a little worse than GS in Tphn but in general much better 
than in the case of the remaining scaling variables (1.5) ~ (1-7). 

Indeed, for the running coupling constant rcl case (1.5) we see in Fig. 5.b monotonous 
fall of (//min) and ((/imin)) with rcj-ef and Xcut respectively, although the large errors at small 
x's allow for a constant fit up to Xref, Xcut — 0.008 yielding 

MO = ((Mmm(0.008))) = 1.677 ± 0.014. (3.7) 

The situation is similar for running coupling rc2 case (1.6) where the constant fit is possible 
up to Xref, Xcut — 0.02 (sce Fig. 5.c) giving 

= ((z.jnin(0.02))) = 2.909 ± 0.025. (3.8) 

In this case, however, one should bear in mind that more "differential" measure of GS - 
X^('^min) - shown in Fig. 3.c does not support hypothesis of GS above x ~ 10~^. 

Finally, in the case of diffusive scaling (1.7) we can hardly conclude that GS is really 
seen; although it is possible to find constant behavior of (Kmin) and {{nmm)) below x ~ 10^^ 
with 

KG = ((Kmin(0.0013))) = 0.449 ± 0.012. (3.9) 

Note, that the errors in Eqs. (3.6) ~ (3.9) are purely statistical (for discussion of 
systematic uncertainties see [3]). 



4 Summary and Conclusions 

In this paper we have applied the method developed in Refs.[3, 4] to assess the quality 
of geometrical scaling of e'^p DIS data on F2 as provided by the combined HI and ZEUS 
analysis of Ref. [24]. In a sense our analysis is in a spirit of previous works [21, 22] 
and especially Ref. [23] where the same set of data has been analyzed by means of so 



- 10 - 



called quality factor. Although the authors of Ref. [23] applied kinematical cuts 4 < < 
ISOGeV^, X < 0.01 our results for scaling parameters given in Eqs. (3.1) and (3.7) - (3.9) are 
in good agreement with their findings. For completeness let us quote their results (note that 
they did not consider logarithmic dependence of r ): ^0 = 1-61 (rcl), 1^0 = 2.76 (rc2) 
and kq = 0.31 (ds). Difference in kq can be explained by applied kinematical cuts, indeed, 
if we take maximal Xcut = 0.01 we obtain ((Kmin(O.Ol))) = 0.301 ±0.006 in agreement with 
[23]. 

Despite the fact that we have been able to find some corners of phase space where 
geometrical scaling in variables (1.5) - (1-7) could be seen, it is absolutely clear that the 
best scaling variable is given by (1.4) (or even by a constant A of Eq. (3.1)), whereas 
diffusive scaling hypothesis is certainly ruled out. This is quite well illustrated in Fig. 4 
where effective exponent Aeg for scaling variable (1.7) changes sign for small Q^. This is 
the reason why in Ref. [23] a cut on low has been applied. Similar argument applies for 
the running coupling rc2 case (1.6) which blows up for small Q^. Because of that x^. ^^^^ 
functions have no minima for very low xi and Xref (points with small x have also small 
Q"^). Therefore the only candidate for scaling variable is running coupling rcl case (1.5). 
Nevertheless, comparing Fig. 5.b with Fig. 6 where we plot results for GS scaling with 
constant exponent A, we see that both by quality and applicability range, the original form 
of scaling variable does much better job than (1.5). Although our results for best values 
of parameters entering definitions of scaling variables (1.5) - (1.7) are in agreement with 
Refs. [21-23] we do not confirm their conclusion that only diffusive scaling is ruled out 
while for other forms of scaling variable geometrical scaling is of similar quality. It is of 
course perfectly possible that the HERA data are not " enough asymptotic" and geometrical 
scaling in one of the variables defined in Eqs. (1.5) - (1-7) will show up at higher energies 
and lower Bjorken x's. 
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